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VARIABLE  COLLISION  FREQUENCY  EFFECTS  ON 
HOSE  AND  SAUSAGE  INSTABILITIES  IN 
RELATIVISTIC  ELECTRON  BEAMS 


I .  INTRODUCTION 

The  resistive  hose  instability  appears  to  be  the  dominant  process 

leading  to  the  observed  disruption  of  propagation  seen  in  many  relativistic 

electron  beam  experiments.*-^  In  these  experiments,  the  beam  is  injected 

into  initially  neutral  gas,  but  beam  ionization  of  the  gas  leads  to  a  rapid 

rise  in  conductivity  a  as  the  beam  passes.  A  number  of  theoretical^  ^  and 
7-13 

computational  models  have  been  developed  to  study  'this  instability,  and 
the  importance  of  the  coupling  between  conductivity  generation  and  beam 
dynamics  has  been  realized  for  some  time.  If  the  plasma  is  weakly  ionized, 


as  is  frequently  the  case,  the  plasma  electron  mobility  is  controlled  by 
the  collision  frequency  vm  for  momentum  transfer  between  plasma  electrons 
and  neutral  molecules.  If  \>m  is  constant,  a  is  proportional  to  the  plasma 
electron  density  ne,  independent  of  the  plasma  electron  temperature  Tg,  and 
usually  increases  monotonically  with  distance  C  5  ct-z  from  the  beam  head. 
Lampe  et  al.^  have  shown  that  analytical  dispersion  relation  estimates  of 
the  hose  growth  rate  can  be  made  if  it  is  assumed  that  the  beam  and  plasma 
current  densities  and  J^,  conductivity  generation  rate  da/dC,  and  beam 
relativistic  factor  y  are  constant  in  C,  and  z.  The  perturbed  conductivity 


cr^ ,  tends  to  follow  the  perturbed  beam  current  J^,  thus  reducing  the 
destabilizing  repulsion  in  the  "normal"  case  where  and  flow  in 
opposite  directions. 
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function  of  E/p  alone.  Furthermore,  p  is  a  constant  for  the  cases 

considered  in  this  paper,  so  v  can  be  regarded  as  a  function  of  E.  In 

most  cases,  Tg  and  E  reach  a  peak  at  the  "pinch  point"  in  the  expanded  beam 

head  and  decrease  with  C  in  the  well-pinched  body  and  tail  of  the  beam.^  ^ 

2  16—18 

Since  9v  /3T  >  0  for  most  weakly  ionized  gases  and  a  =  n  e  /mv  ,  this 

me  e  m 

causes  conductivity  to  increase  more  rapidly  with  C,  after  the  pinch  point 
than  does  ng.  Instability  growth  in  C  is  characterized  by  the  dipole  decay 


length 
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Cl  =  cx1  = 


rttJQ(  -,)a 
2c 


(1) 


where  a  is  the  beam  scale  radius  and  a  s  a  (r  =  0).  The  combination  of 

o  o 

small  n  and  large  v  in  the  front  of  the  beam  can  lead  to  very  rapid  hose 
e  m 

instability  growth  in  this  region.  However,  the  much  larger  ng  and  smaller 

v  in  the  beam  tail  can  cause  hose  growth  rates  to  become  quite  small.  As 
m 

a  result,  the  saturated  transverse  beam  displacement  Y  ( t)  can  approach  a 

Sal 

constant  value  in  the  beam  tail,  contributing  to  the  "plateau"  behavior 
observed  in  many  hose  simulations ^  In  contrast,  Y  (C)  displays 

S3  l 


power  law  growth  if  da/dC  and  v  are  assumed  to  be  constant 
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This  paper  will  also  discuss  a  destabilizing  mechanism  associated  with 

the  perturbed  conductivity^  a.  when  3v  /3E  >  0  or  3v  /3T  >  0.  As  pointed 

out  previously,  reduces  instability  growth  when  is  constant  because 

the  local  conductivity  and  electron  density  increase  in  regions  where  the 

distortion  associated  with  the  instability  increases  the  beam  density. 

However,  the  electric  field  | E^  |  also  increases  in  such  regions,  thus 

increasing  and  decreasing  a.  This  results,  in  some  cases,  in  a 

substantial  increase  in  hose  growth.^  If  v  is  assumed  to  be  a  function 

m 

of  E/p  only,  we  shall  show  that  the  monopole  (unperturbed)  conductivity  in 
the  perturbed  Ampere's  law  can  be  replaced  by 


2 


*  F  3VJ 

•  -  l-rr  ■  ”„<l  -  5).  (2) 

v.  m  y 

Since  q  can  exceed  0.5  for  air,  hose  growth  rates  can  be  substantially 

increased.  This  destabilizing  mechanism  competes  with  the  stabilizing 

effect  of  perturbed  conductivity  described  by  Lampe  et  al.7  Equation  (2) 

represents  one  of  the  major  results  of  this  study. 

In  Sec.  II  of  this  paper  we  derive  equations  for  the  perturbed 

conductivity  and  fields  when  vm  is  a  function  of  E.  Specializing  further, 

we  obtain  analytic  solutions  for  beam  equilibrium  and  hose  instability  for 

the  case  with  constant  dn  /dC  and  v  proportional  to  E0*  (with  q  constant). 

e  m  r 

Hose  simulation  results  obtained  with  the  VIPER  multi-component  simulation 
model7  are  described  in  Sec.  III.  A  sausage  instability  model  using  the 
assumptions  in  Sec.  II  is  presented  in  Sec.  IV,  showing  that  the 
destabilizing  effects  of  the  variable  collision  frequency  are  not  strong 
enough  except  in  extreme  cases  to  give  significant  instability  growth.  The 
implications  of  this  study  are  discussed  in  the  final  section. 
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II.  THEORETICAL  CONSIDERATIONS 

A.  Introduction 

In  this  section,  we  will  derive  some  simple  models  which 

illustrate  some  of  the  important  effects  of  variable  \>m  conductivity  on 

hose  instability.  Much  of  this  work  is  a  direct  extension  of  the 

calculations  of  Lampe  et  al.^  with  constant  v  .  The  spread  mass  beam 

m  r 

dynamics  model  will  be  employed  since  that  model  leads  to  tractable 
dispersion  relations  in  some  situations. 

B .  Conductivity  and  Momentum  Transfer  Frequency 

We  consider  the  case  in  which  the  beam  propagates  in  a  resistive 
medium  which  has  a  scalar  conductivity  given  by 


i 


o  = 


nee 


mv 


(3) 


Many  physical  processes  contribute  to  the  determination  of  a(x,t),  and  a 
variety  of  theoretical  models,  both  simple  and  complex,  have  been  used. 
The  electron  density  ng  is  determined  by  beam-collisional  ionization, 
avalanche  ionization,  recombination,  and  attachment.  However,  we  shall 
concentrate  on  the  "beam  body"  region  where  collisional  ionization 
dominates,  and  ng  is  given  by 


dn 

d^  =  pJb’ 


(A) 


where  u  is  the  production  rate  due  to  beam  impact  ionization  and  is  the 

beam  current  density.  Theoretical  models  are  available  to  prescribe  the 

collision  frequency  \tm  as  a  function  of  T  .  If,  for  example,  the  electron 

distribution  function  f(v)  is  approximately  isotropic  with  a  temperature 

T  ,  then  v  is  given  by 
e  m  & 


_1 _ 4n 

v  (T  )  "  "  3n 
m  e  e 


v  9f 

Ever  ,M.  9v 
mi  1 
1 


dv 


(5) 


where  c  .(E)  is  the  momentum  transfer  cross  section  for  electron  collisions 

19 

with  a  species  i  with  number  density  M...  For  a  Maxwellian  f(v),  Slinker 

and  Ali^  have  compiled  tables  of  \>m(Tg)  for  ^  and  0^  based  on  Eq.  (5). 

Figure  1  shows  \>m(T  ).  Alternatively,  vm  may  be  deduced  from  drift 

20  21 

velocity  measurements  in  discharges,  '  in  which  case,  v  is  known  as  a 

function  of  the  electric  field/gas  density  ratio  E/p.  Over  a  wide  range  of 

validity  Te  may  be  expressed  as  a  function  of  E/p,  so  that  may  in  any 

case  be  regarded  as  a  function  of  E/p.  Electron  heating  is  usually 

dominated  by  ohmic  heating  and  direct  deposition  by  the  beam,  while 

collisions  with  neutrals  result  in  electron  cooling  through  vibrational  and 

21 

electronic  excitations.  Thus 


2 
3  C 


aE 


n 

ev  e 


n  R  (T  ), 
3  n  nv  e' 


(6) 
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where  Cgv  =  1.6x10  ergs/eV,  Sg  is  the  production  rate  due  to  beam  impact 
ionization,  es  is  the  average  energy  of  a  secondary  electron  (=  7.55  eV  for 
air),  and  R  (T  )  is  the  cooling  rate  due  to  collisions  with  neutrals  with 
density  nn<  If  9Te/3t  =  0,  and  ohmic  heating  balances  excitation  cooling, 
then  Te  is  a  function  of  E/p  only.  [This  follows  from  the  fact  that 
a  =  (ng/ p) x( func t ion  of  T  )  for  weakly  ionized  gas.)  Thus  vm  can  be 
expressed  as  a  function  of  E/p. 
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C.  Linearized  Conductivity  and  Field  Equations 

The  hose  instability  is  usually  treated  analytically  by 
partitioning  J^,  a,  Ez>  etc.,  into  monopole  (m  =  0)  and  dipole  (m  =  1) 
quantities, and  linearizing  in  the  dipole  quantities,  which  are  treated 
as  small  perturbations  about  axi-symmetry .  For  example, 

p(r,C,z,0)  =  ffQ(r,C,Z)  +  cr^r,  C,z)cos0.  (7) 

We  will  use  the  0  and  1  subscripts  to  denote  unperturbed  (m  =  0)  and 
perturbed  (m  =  1)  quantities.  In  the  absence  of  avalanche  and 
recombination  processes  the  electron  density  is  given  by 


dnel  , 

dt  _  uJt>r 


If  v)  is  a  function  of  E/p  and  p  is  constant,  then 
m 


n  e 
eo 

ffo  "  m\>  (E) 
m 


2  2 
n  ,e  n  e^  (9v 

el  _  eo  _ m  „ 

CT1  “  mv  (E)  2  9E 

m  mv  t  o; 

m 


fn  .  E,  3m  ) 
el  _  _1 _ m 

7o  n  v  9E  - 
1  eo  m  o) 


The  paraxial  beam  approximation  is  made,  and  the  conductivity  is  assumed 


to  be  sufficiently  large  to  provide  space  charge  neutrality,  so  that 
Er  =  0.  Then  becomes  the  only  significant  component  of  the  vector 
potential , 
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3A 
_ o 
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(10a) 


E  2  E 
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E-  =  E  1 
1  z.1 


3A1 

9T  ’ 


(10b) 


and  Maxwell's  equations  reduce  to  Ampere's  law, 


4ita 


(1  3  3  "  o  3  ) .  4it  .  ni  , 
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f3  1  3  4  Iter  3  4n  4n  a,  3A 
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(lib) 


Using  (9b)  for  and  (10)  for  Eq  and  E^,  Eq.  (lib)  can  be 
rewritten  as 


3  13  4  It  c 

f  E  3m  ^ 

3  ‘ 

4n 

r  n 
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t  c  el 
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Equation  (12),  supplemented  by  Eq.  (8b),  specifies  the  perturbed  potential 

A^  in  terms  of  the  perturbed  beam  current  Jbr  11  is  formally  similar  to 

Eq.  (9)  of  Ref.  7,  which  was  derived  for  the  case  of  constant  v  .  The 

m 

effects  of  beam-generated  perturbed  conductivity  a,  appear  in  two  places  — 
the  last  term  on  the  right-hand  side  (the  "density  term")  and  the  last  term 
on  the  left-hand  side  (the  "field  term").  The  density  term  has  been 
discussed  previously.^  For  the  usual  case  in  which  the  plasma  current 
is  a  return  current  flowing  opposite  to  J^,  the  density  term  is 
stabilizing;  it  leads  to  conductivity  perturbations  which  follow  the  beam 
distortions  and  thus  tend  to  reduce  the  spatial  separation  between  beam  and 
plasma  currents.  The  "field  term",  which  is  discussed  for  the  first  time 
in  this  paper,  does  not  change  the  mathematical  structure  of  the  equation, 
but  in  essence  replaces  the  conductivity  oq( C)  on  the  left-hand  side  by  an 
effective  conductivity 


(13) 


For  "normal"  weakly  ionized  gases  such  as  air,  0  <  q  <  1.  Thus  the 
—  ★ 

(1-q)  term  reduces  and  increases  the  instability  growth  rate.  One  must 

also  note,  however,  that  the  dependence  of  vm  on  Eq  also  directly  modifies 

aQ  in  Eqs.  (12)  and  (13).  In  propagating  beams,  Eq  steadily  decreases,  as 

a  function  of  C,  from  a  maximum  at  the  "pinch  point"  just  behind  the  beam 
22 

head.  Thus,  in  "normal"  gases  v  is  a  decreasing  function  of  £,  which 

tends  to  augment  the  increase  of  as  a  function  of  C.  This  term  thus 
acts  to  speed  up  hose  growth  near  the  front  of  the  beam  and  slow  it  down 
further  back  in  the  beam. 

In  the  next  sections  we  shall  derive  explicit  solutions  of  Eqs. 
(8),  (12),  and  the  equations  for  beam  dynamics  that  illustrate  these 
effects . 

D.  Analytic  Solutions  for  Beam  Equilibrium 

A  number  of  previous  works^  ^  have  derived  analytic  solutions  for 
hose  growth  by  introducing  some  additional  simplifications  and  assumptions 
in  the  assumed  equilibrium  and  the  beam  dynamics.  This  has  proven  quite 
useful,  particularly  for  developing  insight  into  instability  growth  and 
scaling.  In  this  section,  we  shall  calculate  beam  equilibria,  adopting  the 
familiar  assumptions  used  in  these  earlier  works,  but  concentrating  on  the 


new  effects  introduced  by  variable  v  .  In  the  next  section  we  shall  use 

J  m 

the  same  assumptions  to  calculate  hose  growth  analytically. 
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We  assume  that  the  beam  current  1^  is  constant,  and  that  the  beam 


equilibrium  current  density  has  a  Bennett  radial  profile, 


Jbo<r>  = 


.I./na 

b 


bo 


2  2  2 
(1+r  /aV 


2  2  2’ 
(1+r  /a  ) 


(15) 
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where  the  Bennett  radius  a  is  constant.  Lee  has  shown  that,  under 

certain  conditions,  beams  evolve  toward  the  Bennett  profile  as  a  result  of 

23  24 

scattering,  and  this  has  also  been  observed  experimentally.  ’  The 
assumption  of  constant  Bennett  radius  is  reasonable  in  the  beam  body.  It 
follows  from  Eq.  (8a)  that  ngo  also  has  a  Bennett  profile  in  the  beam  body. 
An  approximate  solution  of  Eq.  (11a)  is  then  given  by 


A 
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(r,C) 


2,2. 


A0(Oln[l±V^), 

vl+IT/a  ' 


(16a) 


where  Aq(C)  satisfies 
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(16b) 


and  R  is  a  radius  (usually  R  >>  a)^  ^  where  A  goes  to  zero.  J,  and 
o  J  o  o  6  bo 

nQ  =  uJ^qC,  are  the  axis  values,  independent  of  r.  The  approximation  made 
in  going  from  (11a)  and  (16a)  to  (16b)  is  to  neglect  the  weak  logarithmic 
dependence  of  =  -  3Ao/3C  on  r;  this  also  results  in  the  plasma  current 
density  JpQ(r,C)  having  a  Bennett  radial  profile  with  Jpo(0,C)  5  3^(0- 
We  assume  additionally  that 


vm(E)  =  6 | E |q ,  (17) 

where  8  and  q  are  constants.  Thus,  the  correction  factor  q  in  (13)  or  (14) 
is  constant  and  equal  to  q.  For  many  gases,  this  is  a  good  assumption  over 
a  wide  range  of  values  of  E/p  (with  p  constant).  However,  it  must  be 


9 


f 


EV 


remembered  that  if  |E|  becomes  very  large,  avalanche  breakdown  will  occur, 
thus  invalidating  Eqs.  (8)  and  (17).  If,  on  the  other  hand,  |  E  |  becomes 


very  small,  the  temperature  Tg  will  depend  on  processes  other  than  ohmic 


heating,  the  full  Eq.  (6)  must  be  used,  and  \>m  will  become  independent 


of  E. 


Having  made  these  assumptions,  Eq.  (16b)  can  be  rewritten  as 


,  dA  .1-q  _  2 

>(**)  +  ("  dr)  =  ~  Jbo’ 


where 


2-2  2 
Jia  uJu  e  r  .  R  .,1-q 


where  we  have  also  used  the  fact  that 


if  >  0.  For  constant  JbQ  and  Aq ( Co )  =  0,  Eq.  (18)  has  the  solution 


2  ,t  ,1/k 

’b„ 1  -  (r) 


(21a) 


if  q  =  0,  or 


A0(O  = 


2  v  _a_  _  ,  _  _i_  _  _a_ 

ria  -=  1-q  1-q  T  1-q  .  1-q 

c  bo  o 


(21b) 


if  q  t  0  and  q  i-  1.  We  are  assuming  C  >  C 


We  may  now  consider  the  general  nature  of  this  equilibrium  for 


various  ranges  of  q. 


.*•  ^  ^  . 


bowwow 


Case  ( 1 ) :  q  <  0.  A  singularity  exists  at  some  finite  value  C* 

2 

of  C»  where  A  -  ( na  3,  /c),  E  -»  0,  a  -»  0,  and  3  -»  0.  For 

o  bo  o  o  po 

C  >  t*.  the  only  solution  is  E  -a  =  3  =0.  This  singularity  occurs 

o  o  po  °  J 

because  Eqs.  (3)  and  (1/)  specify  that  a  goes  to  zero  for  E  -»  0.  It  is 

unphysical,  of  course;  Eq.  (17)  loses  its  validity  when  |e|  becomes  small. 

However,  it  can  be  shown  that  the  plasma  current  density  J  decreases  from 

r  po 

its  initial  value  on  a  C  scale  that  is  slow  compared  to  the  scale  length 

2  2 

for  hose  growth  essentially  because  of  the  large  quantity  ln(l  +  Rq  /a  ) 

in  the  definition  of  K,  Eq.  (19).  Thus,  there  is  a  range  of  2,  where  it  can 

be  assumed  that  J  is  constant  which  will  be  seen  to  simplify  the 
po  v  i 

instability  calculation. 

Case  (ii):  q  =  0.  This  case  is  well  known.  Here  Eq  ~  K 

-1/K 

<tq  ~  C,  and  Jp  ~  2.  .  The  plasma  current  decays  to  zero  as  2,  -»  ”,  but 

again,  this  decay  is  slow  compared  to  the  C  scale  for  hose  growth. 

Case  ( iii ) :  0  <  q  <  1 .  Here  Eq  ~  C  aQ  ~  and 

-»  const  as  2  ->  ®.  This  scaling  eventually  breaks  down  at  large  2: 

either  |e|  becomes  so  small  that  Eq.  (17)  becomes  invalid,  or 

alternatively,  recombination  may  become  important  in  limiting  neQ,  in  which 

case  Eq.  (8a)  becomes  invalid.  In  either  case,  JpQ  eventually  decays  to 

zero  asymptotically  at  2,  ->  ®.  However,  it  is  reasonable  once  again  to 

assume  that  JpQ  is  constant  over  the  2,  scale  for  hose  growth. 

Case  (iv):  1  <  q.  This  case  is  very  peculiar,  since  the  effective 

conductivity  ct*  is  negative.  If  a  current  J^o  is  flowing  in  the  gas,  small 

local  increases  in  |  E  j  decrease  the  current!  Gases  do  exist  which  have 

2  6 

this  property  for  certain  ranges  of  E/p.  Christophorou  gives  several 
examples,  such  as  argon  doped  with  C,,^.  Current  flow  is  unstable  in  such 
a  gas,  even  in  the  absence  of  an  electron  beam,  since  the  left-hand  side  of 


Eq .  (12)  is  a  diffusion  equation  with  a  negative  diffusion  coefficient.  It 


is  thus  clear  that  one  cannot  assume  simple  equilibria  of  the  type 

considered  here.  We  shall  not  consider  this  case  further. 

E.  Analytic  Solution  for  Hose  Growth 

1 .  Conductivity  and  Field  Equations 

As  ve  have  seen  in  the  previous  subsection,  if  q  <  1  it  is 

qualitatively  reasonable  to  assume  that  3  is  constant  over  the  C  scale  of 

po 

primary  interest  for  hose  growth.  (However,  we  must  keep  in  mind  that 

eventually  JpQ  does  decay  to  zero  as  £  -»  ®. )  Making  this  assumption,  in 

addition  to  those  discussed  previously,  significantly  simplifies  the  hose 

growth  problem.  We  shall  define  a  current  neutralization  fraction 

f  *  3  /3,  which  is  a  constant.  Equation  (8a)  reduces  to 

po  bo 


n 


eo 


uJboC’ 


(22) 


and  the  density  term  in  the  perturbed  conductivity  in  (12)  is 


a  E 
o  o  n 


el 


eo 


fJbonel 

uJboC 


fn 


el 


pt 


(23) 


Following  Ref.  7,  we  take 


3n 


el 


3C 


-  ('  * c  It) 


el  . 

—  =  yJbl. 


(24) 


Thus,  Ampere's  Law  (12)  takes  the  form 


t1  * c  Is) 


3  13 

3r  r  3r  rAl 


4n 


(,  )  9Ai 

’o  M  8 r 


4n  f.  7  3  f)  T 

c  3C  )  bl 


(25) 


Equation  (25)  is  identical  to  (13)  in  Ref.  7  if  one  replaces 


o  (1-q)  =  a 
o  M  o 
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Spread  Mass  Hose  Instability  Model 


In  order  to  estimate  hose  instability  amplitudes,  the  field 
equation  (25)  must  be  combined  with  a  beam  dynamics  model  which  specifies 
in  terms  of  A^.  Following  Ref.  7  we  assume  that  perturbed  quantities 
have  a  z-dependence  exp(-iSz)  and  employ  the  spread  mass  model  of  Lee."*  In 
this  model,  the  perturbed  current  and  vector  potential  are  assumed  to  be 


rigid  displacements  of  the  equilibrium  profiles:  J. 


Y(dJ^Q/dr)  and 


Aj  =  -  D(dAQ/dr).  The  displacements  Y  and  D  are  related  by 

Y(C)  =  D(C)  (l  +  G(s)),  (26) 

2  2  1/2 
where  s  s  52  /G^,  sgo  s  1 2 n(  1  +  f  )eJ^o(0)/(  ymc)  ]  is  the  on-axis  betatron 

frequency,  and  G(s)  is  Lee's  spread  mass  weighting  function"*  which  for  a 

Bennett  profile  is 


G(s)  =  6s -  s  +  (s  -  s2)[,i  +  lnp^-jjj. 


Inserting  these  relationships  into  (25)  along  with  the  equilibrium 
and  perturbed  radial  dependences  for  J^,  a,  and  A  results  in  the  following 
ordinary  differential  equation  for  the  field  displacement  D(C) 


l1  *  <  ('  *  ic)D  ■  ('  *  *  ♦  <•  y®.  (28) 


where  the  constant  A  is 

q 


with  g  =  e  p/mg. 


We  assume  that  v  -  E  ,  so  that 


s  S 


(30) 


a 

o 


X 

q 


Equation  (28)  reduces  to  the  result  in  Re£.  7  for  q  =  0;  in  that 
case,  a  dispersion  relation  can  be  obtained  by  replacing  (  by  In  (  as  the 
independent  variable.  In  the  present  case,  it  is  not  possible  to  derive  a 
dispersion  relation  in  closed  form,  but  the  hose  growth  resulting  from 
Eq.  (28)  can  be  expressed  in  terms  of  tabulated  functions  or  calculated  by 
numerical  integration. 

3 .  Kummer  Function  Solutions 

The  solution  to  Eq.  (28)  can  be  expressed  purely  analytically  in 
terms  of  confluent  hypergeometric  functions.  Although  not  very  useful  for 
obtaining  numerical  results,  this  method  can  illuminate  some  of  the 
analytical  properties  of  the  instability  for  large  C  in  different  ranges 
of  q. 

25 

Equation  (28)  can  be  rewritten  as  Kummer' s  equation, 


(31) 

(32) 

(33) 

(34) 


The  solutions  are  the  Rummer's  functions  M(a,b,v)  and  U(a,b,w)  with 
arguments  defined  by  (32),  (33)  and  (34).  The  complete  solution  is 


D(C(v))  =  AM(a,b,v)  +  BU(a,b,v),  (35) 

with  A  and  B  arbitrary  constants. 

For  0  <  q  <  1,  which  is  the  range  of  interest  for  weakly  ionized 
air,  the  limit  C  ■*  “  corresponds  to  w  -»  0.  In  this  limit,  we  have 


D(w=0)  =  A  +  B 


rd-b) 

r(l+a-b) 


constant . 


(36) 


This  implies  that  the  hose  amplitudes  in  the  beam  tail  approach  a  finite 

asymptotic  limit  as  £-»«■>,  in  contrast  to  the  case  of  q  =  0,  where  it  was 

0). 

found  in  Ref.  7  that  D(C)  ~  C  and  is  unbounded.  In  principle,  it  is 
possible  to  give  an  analytical  estimate  of  this  bound  by  first  assuming 
that  at  w  =  0,  D  =  1  and  dD/dw  =  0  to  get  A  and  B.  One  can  then  assume 
that  growth  begins  at  a  point  CQ  which  is  sufficiently  small  that 
w(CQ)  a  wq  >>  1.  Then  using  the  asymptotic  expansions  for  M(a,b,wQ)  and 
U(a,b,vQ)  gives 


D(wq)  = 


I  r(b) 

ir(b-a) 


+  i  Jia 
e_  w 


-a 


Uhl 

T(a) 


a-b 


B  w 


(37) 


The  ratio  of  hose  displacements  is  then  given  by 


D(  C  -*  m)  _  D(w  ->  0) 
°(C0)  "  D(vo) 


(38) 


Since  G(s)  is  complex,  w  is  in  general  complex,  so  numerical  solutions  to 
(36-38)  are  not  easy  to  obtain.  It  is  easier  to  simply  integrate  (28) 
numerically,  as  discussed  in  the  next  subsection.  Even  though  the  hose 
amplitudes  have  a  finite  limit,  in  most  cases  of  interest  D(t  ■*  00 )  >> 


Solution  for  q  in  other  ranges  can  also  be  obtained.  For  q  =  0, 


the  transformation  is  invalid,  but  this  case  has  already  been  worked  out  in 


Ref.  7.  For  q  =  1,  the  effective  conductivity  o  ■*  0,  and  the  instability 


is  absolute.  The  case  q  <  0  corresponds  to  having  vm  decrease  with  E,  as 


is  the  case  in  the  Spitzer  conductivity  regime.  In  this  case,  b  ranges 


from  +°°  (for  q  -»  0  )  to  1  for  q  ->  Then  C  -»  ®  corresponds  to  w  ■»  *,  and 


the  solutions  are  (Ref.  22,  Eqs.  13.5.1  -  13.5.2): 


»<->  ■  *  {ft&7  *iina"~a  *  RS  eV‘b) 


The  second  term  is  unbounded  if  Re[G(l+f )/(f-G) ]  <  0,  resulting  in 


unlimited  hose  growth.  This  condition  is  usually  satisfied. 


4.  Direct  Numerical  Integration 


Since  the  purely  analytical  solutions  described  above  are 


cumbersome,  it  is  perhaps  more  useful  to  integrate  the  field  displacement 


equation  (28)  numerically.  This  requires  an  appropriate  choice  for  the 


starting  point  CQ  and  the  constant  X^ .  For  q  =  0,  Lampe  et  al.  used 


Y-o 5  xo  =  -044  V1^’ 


and  estimated  the  return  current  fraction 


‘‘Vf 


with  Dj  =  1 . 


In  order  to  compare  results  with  different  values  of  q,  we  assume 


a  value  of  Xq  and  require  that  at  a  reference  point  C^,  the  dipole  decay 
length  is  equal  to  the  nominal  q  =  0  value,  XQC^.  Combining  (1),  (14),  and 


(30)  results  in 


.VM 


O  \  ■ 


A  v*  V.  J 
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(42) 


CTi  = 


1_ 

1-q 


1-q 


Thus,  the  value  of  X^  which  gives  cx^(C^)  =  is 

_  _L_ 

-  d-q)(X0(„K,  1-t|  .  <i-q)X0CN1"q-  «3) 

The  starting  point  CQ  for  the  perturbation  is  arbitrarily  chosen 

to  be  the  point  at  which  ct«(C  )  =  a  ,  the  nominal  beam  radius.  Given  X 

loo  q 

from  (43),  CQ  can  be  calculated  from  (42)  for  any  chosen  value  of  ct^^). 

Finally,  initial  values  for  Dq  a  D(CQ)  and  dD(CQ)/dC  must  be  chosen;  we 

have  assumed  Re  D  =  Im  D  =  /2/2  and  dD(C  )/dC  =  0. 

o  o  o 

An  example  of  hose  instability  growth  for  various  values  of  q  is 
shown  in  Fig.  2.  The  following  parameter  values  were  chosen:  =  1, 

Xo  =  lf  ^  =  20  a0’  f  =  -°*5>  Sr/Seo  =  °‘5,  and  cT1(^0)  =  a0-  For  q  =  0, 
pure  power  law  growth  is  observed  for  (  >  10  ao  with  a  slope  equal  to  the 
value  <ak  =  1.80  predicted  by  the  dispersion  relation  of  Ref.  7.  For 
q  =  0.50,  hose  growth  is  very  rapid  in  the  front  of  the  beam  but  flattens 
out  into  a  plateau  for  C  >  100  aQ.  This  qualitative  behavior  is  seen  for  a 
wide  range  of  Xq  values  and  becomes  more  pronounced  as  q  approaches  1. 
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III.  HOSE  INSTABILITY  SIMULATIONS  VITH  VARIABLE  V  MODELS 

m 

A.  Introduction 

The  hose  growth  behavior  shown  in  Fig.  2  has  been  observed  in  a 

variety  of  simulation  models  including  VIPER,  SIMMl,^  SARLAC,^ 

9  11 

EMPULSE,  and  PHLAP.  These  codes  allow  J^,  J  ,  and  ng  to  vary  with  r,  £,, 

and  z  and  permit  J^,  J^,  and  a  to  have  different  radial  profiles.  In  this 

section,  results  from  the  VIPER  multi-component  simulation  model  will  be 

presented  for  cases  with  \>m  =  vm(E/p).  A  description  of  the  VIPER  model 

and  some  examples  with  a  fixed  v  conductivity  model  are  contained  in 

m 

Ref.  7. 

B.  VIPER  Simulation  With  v  ~  Eq 
-  m 

The  case  described  in  Sec.  II  can  be  simulated  by  using  a  momentum 
transfer  frequency  of  the  form 

vm(E/p)  =  ♦  vJ(E/p)q.  (44) 

The  floor  on  is  used  to  avoid  unrealistically  large  conductivity  in  the 

weak  field  regions  of  the  beam.  Avalanche,  recombination,  and  beam 

expansion  due  to  scattering  have  been  turned  off  in  the  simulation,  but 

and  do  vary  with  z  and  C.  The  principal  quantity  of  interest  is 

Y  (C)>  the  maximum  or  "saturated"  displacement  of  the  beam  slice  at  a 
S3 1 

distance  C  from  the  beam  head  as  it  propagates  forward  in  z. 

Figure  3  plots  Y  (C)  for  a  series  of  VIPER  simulations  in  which 

S3  L 

the  exponent  q  in  (44)  was  varied.  The  initial  perturbation  began  at 


C  =  a  ,  and  \>  was  chosen  to  give  a  reference  momentum  transfer  frequency 
o  o  m  6  m  j 

12  -1 

of  1.4  x  10  sec  at  E/p  =  3  kV/cm-atm.  As  q  is  increased,  Y  , ( C)  shows 

S3  l 


the  qualitative  behavior  seen  in  Fig.  2.  Higher  values  of  q  give  more 
pronounced  plateau  behavior. 


>»*5 


l:5:| 
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C.  VIPER  Simulations  With  More  Elaborate  Conductivity  Model 

18 

A  more  elaborate  conductivity  model,  based  on  rate  equations  and 

v  (T  )  calculations  from  Slinker  and  Ali*7  and  Ali,"^  has  also  been  used 
m  e 

21 

with  VIPER.  The  model  assumes  a  Tg(E/p)  relation  deduced  by  Ali: 


tiVl 

\m 


Vi1 

$ 

& 


•sa 


T  =  T  +  0.1(E/p) 
e  eo 


with  E/p  in  volts/cm-torr ,  T  in  eV,  and  T  =  0.15  eV.  The  total 

e  eo 

collision  frequency  is  the  sum  of  vm  and  the  electron-ion  contribution. 

Figure  4  shows  the  maximum  hose  displacement  (saturation  had  not 

occurred  by  the  time  of  the  graph)  as  a  function  of  C  for  a  10  kA,  5  mm 

beam.  The  three  curves  correspond  to  three  VIPER  runs.  In  A,  the  full 

dipole  conductivity  terms  in  (9b)  are  retained.  In  B,  the  dipole 

conductivity  is  set  equal  to  zero.  In  C,  only  the  contribution  to  a ^  due 

to  ng^  is  retained.  As  seen  in  the  figure,  this  part  of  the  dipole 

conductivity  is  stabilizing. 7  The  inclusion  of  the  3v  /9E  term,  Case  A, 

mo 

is  seen  to  be  destabilizing  in  accordance  with  theory. 


m 

w 


Finally,  Fig.  5  shows  the  on-axis  conductivity,  a,  electric  field 
_  ★  _ 

and  fraction  a  /a  as  a  function  of  C.  Since  most  of  the  hose  growth  occurs 
in  the  region  50  <  C/aQ  <  300  where  a*/ a  =  0.6  is  approximately  constant, 

Eq .  (12)  suggests  that  lnY(Case  C)/lnY(Case  A)  should  be  approximately  0.6. 


I 


Figure  4  confirms  this  expectation. 


IV.  SAUSAGE  INSTABILITY  FOR  v  -  Eq 

HI 

A.  Introduction 

The  resistive  sausage  instability  is  excited  by  the  repulsive 

force  between  the  beam  current  and  return  current.  This  is  an  axi- 

symmetric  (m  =  0)  mode  in  which  the  beam  expands  or  contracts  self- 

similarly.  For  an  infinite  beam  with  constant  conductivity,  the 

instability  is  excited  when  the  ratio  -I  /I.  =  -f  exceeds  0.69  for  a 

p  b 

28 

Bennett  profile.  However,  the  instability  to  our  knowledge  has  never 

been  unambiguously  observed  in  simulations  or  experiments  in  which  beam 

impact  ionization  is  the  dominant  conductivity  generation  mechanism. 

28 

Lampe  and  Joyce  have  shown  that  for  the  assumptions  used  '.n 
Section  III  (f,J^Q,  and  doQ/dC  constant),  and  if  in  addition  q  =  0,  the 
inclusion  of  dipole  conductivity  effects  leads  to  the  following  instability 
threshold  condition  on  the  return  current  fraction  f: 

(2X  .  /3)2 

■  nr  >  -sr  ~i  -  <«> 

o 

The  function  H(a,XQ)  is  a  slowly  varying  function  which  approaches  0.8  when 

29 

the  anharmonic  damping  coefficient,  a,  equals  Lee's  value  of  0.7,  and 

H(a  -*  <E.Xo)  -*  1.  If  one  estimates  f  from  (41),  the  instability  condition 

is  not  reached  unless  >  4,  a  situation  almost  never  observed  in  beam 

28 

simulations.  Thus,  Lampe  and  Joyce  concluded  that  the  sausage 
instability  should  be  stabilized  in  almost  all  conceivable  cases  for  beams 
propagating  in  initially  neutral  gases. 

In  the  next  section,  we  extend  the  analysis  of  Ref.  28  to  the  case 
where  ~  (E/p)q.  This  introduces  the  destabilizing  dipole  conductivity 
term  in  Ampere's  Law  (12)  and  reduces  the  threshold  return  current  fraction 


B.  Derivation  of  Sausage  Dispersion  Relation 

28 

Lampe  and  Joyce  used  an  envelope  equation  approach  to  treat  the 


linearized  sausage  instability.  The  perturbed  beam  and  plasma  currents  SJ 


and  SJ  are  assumed  to  have  the  form 
P 


SJ 


J(t,z) 


i  2/  2 

1-r  /a 

o 


r  i  2/  2,3 
(1+r  /aQ] 


(47a) 


and  the  perturbed  vector  potential  is  given  by 


SA 


D(C,z) 


i  2.2 

1-r  /a 
_ o 

2  2 " 
1  +  r  /a^ 
o 


(47b) 


The  unperturbed  conductivity  is  given  by  (9a),  and  the  perturbed 

conductivity  has  the  radial  dependence  given  in  Eq.  (19)  of  Lampe  and 
28 

Joyce.  The  perturbed  envelope  equation  can  be  written  as 


2  a2jb 

532  _ b 

6  3z2 


(2/3)Jp 

T+f 


(48) 


where  we  have  changed  the  notation  from  that  of  Ref.  25  to  be  consistent 
with  the  rest  of  this  paper.  The  average  betatron  frequency  8^  =  = 

aQ  ^(I^(l+f )/IA)1/^  where  1^  =  ymc'Ve  is  the  Alfven  current. 

Using  the  assumed  radial  dependence  for  the  various  quantities  in 
Ampere's  Law  results  in 


t1  *  C  fc)Jp  ■  fV  (‘  ‘  C  fth^om-c,)  |2.  (49) 


To  close  the  system,  we  note  that 


na 


D  = 


o  (~  ~ 

=  J.  +  J  , 

2c  b  p) 


(50) 
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and  Fourier  analyze  in  z,  assuming  a  z-dependence  exp(-  iSBftz).  This  gives 

3  (1+f) 

Jb  12  *  U  2 

-Q  -  i aO  +  1 2  +  j 

Combining  (49),  (50),  and  (51)  results  in 

* <  k  -  if-  <52> 

We  have 

2  1 
na  3 — 

2^  ffQ(r=0) (1-q)  =  XqC  q,  (53) 

leading  to  the  following  sausage  growth  equation: 

1 

(i  *  <  3C  -  £s(=>)5  ■  -  \ (i*C  fcK‘q  jjf,  (54) 

with  g( Q)  defined  by  (51)  and  Aq  defined  by  (53). 

C.  Numerical  Results  for  Sausage  Instability  Growth 

In  this  section,  numerical  solutions  to  the  sausage  instability 

growth  equation  (54)  are  presented.  The  method  is  similar  to  that 

discussed  in  Sec.  II.  E.  for  the  hose  instability.  Figure  6  plots  the 

sausage  amplitude  D(C)  with  parameters  q^,  XQ,  C^,  CQ,  and  f  the  same  as  in 

Fig.  2.  The  frequency  Q  is  1.25  (Q/Q^  =  0.88)  which  is  close  to  the  value 

®max  ^or  max*mum  growth,  and  the  damping  coefficient  in  (51)  is  Lee's 
29 

suggested  value  a  =  0.7.  ®max  varies  weakly  with  X  but  is  strongly 
dependent  on  a.  Figure  6  shows  no  sausage  instability  even  for  q  =  0.5 
for  this  case  in  which  f  has  its  nominal  value  (f  =  -0.5  for  q^  =  1  and 
qX  =  1). 


■  *<«>V  <51> 
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Lampe  and  Joyce  point  out  that  sausage  instability  can  be 


generated  in  the  q  =  0  limit  if  q^  exceeds  the  limit  in  (46).  Figure  7 
plots  D(C)  for  a  case  in  which  -f  has  been  raised  to  0.85  (q^  =  5.67).  For 
q  =  0,  the  beam  is  very  weakly  unstable,  as  expected  since  q^  is  slightly 
above  the  threshold  given  by  (46).  As  q  is  increased,  significant  sausage 


growth  appears.  The  critical  of  q^  value  for  the  onset  of  hose  instability 

22 


drops  to  q^  =  3.2  (f  =  -0.76)  when  q  =  0.5.  Since  Sharp  and  Lampe  quote 
a  maximum  value  of  q^  of  3  to  4  in  typical  beams,  this  reduction  in  the 
threshold  for  sausage  instability  may  allow  beams  to  be  weakly  unstable  if 
q  is  sufficiently  large.  The  hose  instability  linear  amplitudes  will  be 


substantially  larger,  however.  Although  sausage  tends  to  grow  out  of 

28 


perturbations  which  are  larger  than  those  for  hose,  any  experiment  which 
is  sausage-unstable  is  likely  to  have  even  more  violent  hose  disruptions. 

D.  Axi-Symmetric  Hollowing  Instability 

A  violent  axi-symmetric  hollowing  instability  has  been  observed  in 


f  .  15,27,30,31  ,  , 

many  particle  simulations  of  propagating  beams.  Joyce  and 


15 


Lampe  showed  that  this  instability  is  driven  by  avalanche  ionization  in 


the  high  E/p  region  near  the  front  of  the  beam  and  proposed  that  the 


instability  is  triggered  if  (E/p)^^^  >  130  kV-cm  *-atm  *  s  (E/p)  ^  , 


Ekdahl  et  al.  have  observed  violent  hollowing  instability  on  the  IBEX 
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experiment  for  beams  propagating  at  pressures  of  10-80  torr.  Simulations 
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by  Godfrey  have  reproduced  the  hollowing  instability  at  pressures 


somewhat  above  the  80  torr  threshold  and  have  confirmed  that  the 

15 


instability  threshold  is  close  to  the  Joyce  and  Lampe  prediction. 

Although  no  detailed  theory  of  the  avalanche-driven  hollowing 
instability  has  appeared,  the  instability  is  believed  to  be  due  to  an 


enhancement  of  a(r)  and  J  (r)  near  the  higher  E^  region  around  the  beam 
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axis.  The  analysis  is  complicated  by  the  fact  that  no  Bennett-like 

equilibrium  exists  when  return  current  flow  near  the  axis  is  enhanced,  so 

the  simulations  appear  to  be  attempting  to  achieve  a  new  equilibrium  at  the 

same  time  as  the  instability  is  being  generated. 

Simulation  studies  to  date  have  been  made  almost  exclusively  with 

constant  v  conductivity  models.  The  inclusion  of  a  v  (E/p)  model  could 
m  J  m 

affect  the  instability  in  several  ways.  Since  the  instability  is  triggered 

ef  f 

in  a  high  E  region,  \>m  and  would  undoubtedly  be  smaller  for  a 

\>m(E/p)  conductivity  model,  thus  increasing  E^Q  and  increasing  the  growth 

rate.  On  the  other  hand,  the  radial  profile  of  a(r)  would  be  less  peaked 

on-axis  due  to  the  higher  vm(E/p)  near  the  axis  which  would  tend  to 

increase  the  threshold  for  instability.  It  is  not  clear  which  of  those 

effects  would  dominate.  However,  since  the  instability  is  driven  primarily 

by  the  strong  nonlinearities  in  the  avalanche  rate  coefficient  for 

E/p  =  (E/p)ct_^  t ,  it  is  likely  that  the  variable  effects  will  not 

change  the  threshold  significantly.  Studies  with  the  SIMM0  particle 
15  18 

simulation  using  the  same  v^fE/p)  model  contained  in  the  VIPER  code  do 
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not  indicate  a  significant  change  in  the  hollowing  instability  threshold. 


V. 


CONCLUSIONS 


The  growth  of  the  resistive  hose  instability  in  a  relativistic 

electron  beam  is  strongly  dependent  upon  the  details  of  conductivity 

generation  by  the  beam.  We  have  considered  cases  in  which  beam  impact 

ionization  is  the  dominant  conductivity  generation  mechanism  and  the 

collision  frequency  v  is  a  function  of  E/p,  the  ratio  of  the  electric 

field  to  the  neutral  gas  density.  The  introduction  of  variable 

introduces  two  important  new  features.  First,  since  3\>m/3(E/p)  >  0  for 

normal  weakly  ionized  gases,  and  since  the  electric  field  usually  decreases 

with  distance  C  from  the  beam  head,  the  local  dipole  decay  length 
2 

ct^  =  np(r=0)a  /2c  tends  to  be  shorter  in  the  front  of  the  beam  but  longer 
in  the  beam  tail  compared  with  values  usually  estimated  from  constant 
collision  frequency  models.  As  a  result,  hose  amplitudes  grow  very  rapidly 
near  the  front  of  the  beam  but  taper  off  and  approach  a  constant  value  in 
the  beam  tail.  The  second  feature  is  a  destabilizing  effect  which  is  due 
to  localized  decreases  in  the  perturbed  or  dipole  conductivity  driven  by 
the  increase  in  collision  frequency  resulting  from  perturbations  in  the 
electric  field. 

The  destabilizing  effect  is  illustrated  in  Eq .  (12)  which  shows  that 
the  term  multiplying  3A/3C  in  Ampere's  Law  is  (4naQ/c) ( 1-fE/v^) ( Sv^/ 3E) ) . 
Thus,  the  beam  behaves  as  if  the  monopole  conductivity  a ^  were  replaced  by 
a  (l-(E/v  )(3v  /3E))  =  a  (1-q).  Since  q  >  0.5  in  some  cases,  the 
destabilizing  effect  can  be  quite  strong.  For  cases  in  which  q  is  constant 
and  otherwise  the  assumptions  of  Ref.  7  are  justified,  a  differential 
equation  for  the  hose  displacement  D(  C)  is  derived.  Solutions  to  this 
equation  in  terms  of  confluent  hypergeomet l ic  functions  are  given,  as  are 
numerical  solutions.  Plots  of  D( C)  show  the  pattern  of  rapid  hose  growth 
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followed  by  a  plateau  which  uecomes  more  pronounced  as  q  approaches  1.  As 
q  -*  0,  the  D(C)  plots  approach  the  pure  power  law  behavior  described  in 
Ref.  7. 

Simulations  of  hose  instability  growth  using  the  VIPER^  model  allow  a 

more  self-consistent  treatment.  In  this  model,  the  monopole  and  dipole 

beam  and  net  current  densities,  etc.  can  vary  with  C  and  propagation 

distance  z.  The  two  major  features  of  variable  collision  frequency 

described  above  are  observed  in  the  simulations  and,  in  fact,  were  noticed 

in  simulations  before  the  analytical  models  were  obtained. 

Variations  in  momentum  transfer  frequency  can  also  affect  the  sausage 

28 

instability.  We  have  extended  earlier  calculations  of  Lampe  and  Joyce  to 

include  these  effects  and  find  that  although  variable  \>  dipole 

m 

conductivity  effects  lower  the  threshold  for  instability,  strong  sausage 
instability  growth  is  unlikely  except  in  extreme  cases.  Such  cases  would 
be  violently  hose  unstable  as  well. 
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3  Saturated  hose  amplitude  for  q  =  0,  0.1,  0.2  and  0.3.  Results 
obtained  from  VIPER  simulations. 
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4  VIPER  simulations  showing  the  effects  of  dipole  conductivity  on 
hose  instability  growth.  Case  A  has  all  the  terms  included 
(Eq.  9b).  Case  B  has  a]  =  0.  Case  C  has  only  the  direct 
production  contribution  to  : 


a,  =  <r  — 
1  on 
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